function nz = func_nz_from_nxy(nxy, CONSTS)

    eps_a = CONSTS.eps_a;

    nz = sqrt_from_negative_num(eps_a-nxy.^2); 
    
end